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The p meson leptonic decay constant at the physical pion mass is calculated with chiral fermions. 
Overlap valence quarks are used on the L® x T = 48® x96Y/ = 2 + l domain-wall fermion 
configurations with light sea quarks at the physical point and the spatial extension of 5.5 fm. The 
mass of the vector meson p is determined to be = 775.9 ± 6.0 ± 1.8 MeV, and the valence up 
sigma term is extracted to be = 9.82(6) MeV from Feynmann-Hellman theorem. By adopting 
a new method to reduce the noise of the wall-wall correlation functions of hadrons, we obtain the 
p leptonic decay constant fp = 208.5 ± 5.5 ± 0.9 MeV, which is in very good agreement with the 
experimental value fp = 209.4 ± 1.5 MeV derived from the r decay r pv^. 

PACS numbers: 12.38.Gc, 13.20.Jf, 14.40.Be, 12.39.Fe. 


The standard quantities that are used to test lattice 
QCD against experiments include hadron masses, pseu¬ 
doscalar meson decay constants and //<-, and Bk [11. 
On the other hand, the lattice results are getting pre¬ 
cise enough, such as fD„ Si, so that any significant 
deviation from experiments could be taken as a signal 
for physics beyond the standard model. Now that lattice 
configurations with Nf =2-1-1 flavors are available for 
physical pion mass and large volume, the list of physical 
observables should be expanded to have a broader exam¬ 
ination of lattice calculations against experiments. The 
obvious choices include the vector meson decay constants 
for p, K*, D* and D* and the vector to pseudoscalar ra¬ 
diative decays V -T Pj, etc. In this work we shall report 
our analysis of the leptonic decay constant fp± of the 
charged p meson, which describes the coupling of and 
lepton pair W through the charged weak current interac¬ 
tion. 

Gauge configurations oi Nf = 2-1-1 domain-wall 
fermions with large spatial volume and physical pion 
mass have been generated by the RBC & UKQCD Col¬ 
laborations [I|. This work is based on the 481 gauge 
ensemble with lattice size L® x T = 48® x 96 [l| . The lat¬ 
tice spacing has been determined to be = 1.730(4) 
GeV, such that the spatial extension of the lattice is 
approximately La ^ 5.5 fm. The light sea quark 
mass is set to give the pion mass = 139.2(4) 

MeV. For the valence quarks, we adopt the overlap 
fermion action, which is another realization of the chi¬ 
ral fermions on the lattice. The low-energy constant 
Aniix, which measures the mismatch of the mixed valence 
and sea pion masses between the domain-wall fermion 
and the overlap fermion, is shown to be very small [J. 
Since overlap fermion accommodates the multi-mass al¬ 
gorithm and the eigenvectors are the same for different 
quark masses, we use 1000 pairs of eigenvectors plus the 
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TABLE I: The table lists the pion masses m,r, the pion decay 
constants fn, and the masses of p at different bare valence 
quark masses. 


(val) 
arUg ’ 

0.00170 

0.00240 

0.00300 

0.00455 

0.00600 

0.02030 

mTT (MeV) 

114(2) 

135(2) 

149(2) 

182(2) 

208(2) 

371(1) 

/.(MeV) 

130.3(9) 

131.0(9) 

131.6(8)) 




mp(MeV) 

773(7) 

775(6) 

779(6) 

784(5) 

789(5) 

836(3) 


zero modes for deflation in calculating quark propa ga- 
tors for several masses on 45 configurations (see Ref @ 
for details). The bare mass parameters are chosen as 

= 0.00170,0.00240,0.00300,0.00455,0.00600 and 
0.02030, which give the pion mass ranging from 114 to 
371 MeV. In this way we can discern the chiral behaviors 
of the mass and the leptonic decay width of the p meson. 

We first extract the decay constant of pion according 
to the partially conserved axial current relation 

rnif^ = {rriu + TOd)(0|u75d|7r). (1) 

For overlap fermions, the quark mass renormalization 
constant and the renormalization constant Zp oi the 
psuedoscalar density u'j^d satisfy the relation Z^Zp = 1. 
Thus the calculation of can be carried out straightfor¬ 
wardly by taking rriu = rrid = after the matrix 

element in Eq. (ED is obtained from the {A4P) and (PP) 
correlation functions. The pion masses and decay con¬ 
stants are listed in TableHl Since the pion masses are very 
close to the physical pion mass, we perform simply a lin¬ 
ear interpolation in and get /,r = 131.3(6) MeV at the 
physical pion mass m,r = 139.5 MeV, which agrees with 
RBG&UKQCD’s result = 131.1(3) MeV on the same 
lattice and their final theoretical prediction = 130.2(9) 
MeV [l|. RBC&UKQGD also calculates on a larger 
lattice with a smaller lattice spacing, = 2.359(7) 
GeV, with a result = 130.9(4) MeV. Their /^s on the 
two lattices imply very small finite a artifacts. This com¬ 
parison can be taken as a calibration of our formalism. 
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FIG. 1: The effective mass plateaus of p at m-^ = 208(2) 
MeV. The blue and black points are from the correlation 
functions C{r = 0, t) and C(r = 4.58a, t), respectively. The 
red points are from the mixed correlation function = 

C(0,t) + ujC(4.58a, t) with the mixing parameter cj = 10. 
The red band shows the fitted mass in the time range tja 
[3,13). The blue line and the black line are the energies of 
two non-interacting pions 2i5.n- with the momentnm modes 
pa = 27r/L(0, 0, ±1) and pa = 2njL{Q, ±1, ±1), respectively. 

When calculating the two-point functions in the vec¬ 
tor channel, we first perform Coulomb gauge fixing on 
the 45 gauge configurations. The quark propagators are 
computed using a wall source, thus the source operators 
of the vector can be expressed as 

(^) = “(y’ ^)7id(z, t). (2) 

y.z 

For the sink operators of the vector, we use the spatially 
extended operators Oy_i(x, <;r) by splitting the quark 
and antiquark field operator with different spatial dis¬ 
placements r, namely, Oy_i(x, <;r) = u{x,t)jid{x + r,t). 
Subsequently the two-point functions C{t,r) with differ¬ 
ent spatial separation r are calculated as 

= ^ (0|Ow(x,t;r)o(7f (0)|0), (3) 

x,2,|r|=r 

where W is the number of rs that satisfy |r| = r. The ef¬ 
fective mass plateaus of C(r, t) with r = 0 (blue points) 
and r = \/^a = 4.58a(black points) at = 208(2) 
MeV are plotted in Fig. [TJ It is seen that the plateaus 
are converged to a uniform effective mass in large time 
range. The difference of the plateaus in the short time 
range shows the r-dependence of the contamination from 
higher states. In order to reduce the higher state contam¬ 
ination, we linearly combine the two correlation functions 
as C'mix(i) = C{0,t) + ujC{4:.58a,t) with an optimal mix¬ 
ing parameter w = 10, by which we can get a very flat 
effective mass plateau starting from t/a = 3, as shown 
in the figure (red points). We fit Cniix(t) using a single¬ 
exponential form in the time range t/a € [3,13] and get 


my Cl = 0.456(4) (plotted as a red band), where the error 
is statistical and is obtained through a jackknife analysis. 

Through the large W analysis, it is suggested by 
Jaffe Q that the p meson is an ordinary meson like 
the Feshbach resonance which exists as a confined qq 
state in the continuum with zero width in the large W 
limit. Its Breit-Wigner width manifests the coupling of 
the bound state and the scattering tttt states, which van¬ 
ishes as Nc goes to infinity, while its pole mass is insen¬ 
sitive to the strength of the coupling Q. In the contin¬ 
uum, the resonance parameters of p, nip = 775 MeV and 
Fp Ri 150 MeV, have been extracted in the J = 1 tttt 
scattering channel. On a finite lattice, the eigenstates 
of the QCD Hamiltonian in the p channel are actually 
the superpositions of the discretized interacting tttt scat¬ 
tering states and the would-be qq confined state, and 
the energy shifts of these eigenstates from the interacting 
TTTT states or the qq state encode the resonance proper¬ 
ties. Many lattice calculations [MU have been carried 
out to extract the resonance parameters of p according 
to the Liischer method [l2|. It is known that the lat¬ 
tice spectrum in the p channel should show an avoided 
level crossing behavior when the energy of tttt states ap¬ 
proaches to the mass of the resonance by varying the 
lattice size or the relative momentum of the tttt states, 
and the mixing of the would-be qq confined state and the 
nearby scattering state gives two eigenstates. For the 
spatial extension La ^ 5.5 fm of our lattice and in the 
case of m,r = 208 MeV, the two lowest energy thresh¬ 
olds of the P— wave tttt states are 2i?,r(001) = 614 MeV 
(shown in Fig.[T]as a blue line) and 2£'.n.(011) = 761 MeV 
(shown in Fig. [T] as a black line) with the relative mo¬ 
menta pa = 27r/L(0, 0, ±1) and pa = 27r/L(0, ±1, ±1), 
respectively. Since 2i?^(001) is far from the expected p 
mass, the corresponding tttt state should mix little with 
p and therefore have an energy close to 2i?^(001), but 
we do not observe this state in the vector correlation 
function. The reason for its absence might be due to 
the implementation of the Coulomb wall-source opera¬ 
tor which suppresses the scattering states with non-zero 
relative momenta (we will discuss this in depth in an¬ 
other publication). However, 2£'.n.(011) is very close to 
the expected p mass, so the mixing of the correspond¬ 
ing TTTT state and the qq state can be sizable. According 
to the feature of the avoided level crossing, the mixing 
gives two states with nearly degenerate energies around 
TOp, which, however, cannot be resolved in our analysis. 
In this sense, we expect that the plateau in Fig. [1] illus¬ 
trates the effective energy of the two adjacent states, and 
will take the fitted energy value mya = 0.456(4) as an 
approximation of the p mass in the lattice unit. 

We take a similar analysis procedure for the correlation 
functions at other pion masses, and the extracted masses 
of p are listed in Table U and are also plotted in Fig. [5] 
with respect to m^, by which the chiral behavior of nip 
can be investigated. A new chiral extrapolation formula 
based on a modified MS regulator and a power counting 
scheme [l^ gives nip{niT^) = mp(0) -I- cim^ -I- C 2 m^ -I- 
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FIG. 2: The p mass nip = 775.9 ± 6.0 MeV at the physical 
point is interpolated by nip{m-,r) = nip(O) + cim^. The red 
band shows the error of the interpolation. The experimental 
value of nip = 775 MeV is also plotted for comparison. 


In where the log term is due to self energy. 
However, from Fig. [5] it is seen that nip is very linear 
in for m.^. ranging from 114 MeV to 371 MeV. A 
correlated jackknife analysis using the form 

?71p(TO^) = TOp(O) + ClTO^, (4) 

gives mp(0) = 766(7) MeV and ci = 0.505(3) GeV"^ 
with ^jd.o.f = 0.13. The p mass at the physical 
is TOp = 775.9 ± 6.0 ± 1.8 MeV, where the second er¬ 
ror is due to the 0.23% uncertainty of the lattice spac¬ 
ing. Our data cannot discern higher order terms in m-Tr. 
We would like to point out that our study is carried out 
for the first time in the chiral region around the physi¬ 
cal point and with chiral fermions, althou gh t here have 
been many lattice studies on this topic . We note 

that Cl is precisely determined, and serves potentially as 
a constraint on the chiral perturbation study of p. Fur¬ 
thermore, Cl is exactly the valence or connected insertion 
part of the tt/o sigma term from the Feynman-Hellman 
theorem dmp/dmj = cim^, since the sea is 

fixed in our partially quenched calculation of nip. One 
can determine the disconnected part from a direct cal¬ 
culation of the mipij} matrix element in the disconnected 
three-point correlator. From the fitted ci in Eq. o , we 
find = 9.82(6) MeV. 

The leptonic decay constant fp is also an important 
characteristic quantity of p. Most of previous lattice 
studies on L, were performed at relatively heavy pion 
masses [isl - l^ . Here we would like to extract the decay 
constant of charged p at the physical point. For charged 
p, for example, fp- is defined by 

(0|J^")(0)|p"(p,C)) = mpfp-CpipX), (5) 

where Jp \x) = (u'^pdf^x) is the charged vector current 


FIG. 3: The red points show 4>i(r) which depicts the fall- 
off of C{r, t) when r increases. The parameterization of 
'l>i(r) is plotted by the curve. The blue points are the ra¬ 
tios Iiirf)/I{oo) at different r^. 


and ep{p,f) is the C-th polarization vector of p with 
C=l,2,3. 

The spatial components of Jp ^ (x) are actually the op¬ 
erators Ov,i(x;r = 0); therefore, the matrix element de¬ 
fined in Eq. ([5]) can be extracted from C(r = 0,1). The 
key challenge is to divide out the matrix element of the 
wall source operator (0|o[7i | V(p, C)). Usually this ma¬ 
trix element can be derived by calculating the wall-wall 
correlation function 

cW(i) ^ = 

r i 

= ^J^(0lOvp(x,t;r)ol^^^’U0)l0}, ( 6 ) 

x,r,2 

where the last equation uses the definition of C (r, t) in 
Eq. (j31). However, a very large statistics is required to 
obtain a satisfactory signal-to-noise ratio for this kind 
of correlation function. The reason for noisy 
is analyzed as follows. Using the spectral expression 
C{r, t) =J2 4?„(r)e“^"‘, when 1 —>• oo one has 

i 

^ (7) 

r 

In practice, we calculate C (r, t) for r ranging from 0 to 
10a and observe the profile of 4’i(r) for t = 7a where all 
C{r,t) are almost saturated by the ground state. The 
$i(r) at = 208(2) MeV (normalized as $ 1 ( 0 ) = 1) 
is plotted in Fig. [3] as red points. It is seen that <i>i(r) 
damps rapidly with r and can be parameterized as 

$i(r) = $i(0)e”(^) , (8) 

with the parameters a = 1.60 and rg = 5.88a. The curve 
illustrates this parameterization in the figure. We check 
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TABLE II: The renormalization constant Za obtained at 
different pion masses, which also gives Zy by the relation 
Za = Zy for overlap fermions. The renormalized decay con¬ 
stants fp are also listed in the table. 


m,r(MeV) 

114(2) 

135(2) 

149(2) 

182(2) 

208(2) 

371(1) 

> 

Za 

1.103(4) 

1.103(3) 

1.104(2) 

1.104(2) 

1.105(1) 

1.105(1) 

o 

/p(MeV) 

206(7) 

208(7) 

211(6) 

215(5) 

217(5) 

223(3) 

a. 


this at other pion masses and find $i(r) is similar for 
all the cases and is very insensitive to m^n-. This means 
that when calculating the wall-to-wall correlation func¬ 
tion the C{r,t)s ( see in Eq. ®) with very large 

r contribute only noise and make (t) very noisy. In 
order to circumvent this difficulty, we introduce a cutoff 
Tc to exclude the contribute of C{r,t)s with r > from 
(t) and use the correlation function 

NrCir,t), (9) 

r<ra 



Exp(( 




0 0.03 0.06 0.09 0.12 0.15 

m/(Gev2) 


FIG. 4: The decay constant fp obtained at different pion 
masses. The curve shows the linear interpolation in terms of 
m^. The extracted from the r leptonic decay is also 

plotted as an asterisk. 


to approximate Letting /i(r') = drr^(^i(r), 

one can see that the ratio (t) can be de¬ 

picted by the ratio /i(rc)//i(oo) at large t. The ratio 
hire)/hioo) using the parameterization above is also 
plotted in Fig. [31 It approaches to 1 beyond re = I5a 
and is equal to 0.995 at Vc = 20a, whose deviation from 
one is already much smaller than the statistical error. So 
we take (20a, t) as a satisfactory approximation of 
( 7 (w) throughout this work. 

With this prescription, we carry out the joint fit to the 
correlation functions 

C(0,t) = 

n 

C'("')(20a,t) 

n 

where is the matrix element of the wall source op¬ 
erator between the vacuum and the n-th state, and fn 
is the bare decay constant of the n-th state according to 
the definition in Eq. ([5]). In practice, two exponentials 
are used in the fit and fi is taken as the bare decay con¬ 
stant of p (the second term is introduced to account for 
the contamination of higher states). Since /i is sensitive 
to the value of rni, we adopt the single-elimination jack¬ 
knife analysis procedure as follows. On each jackknife 
re-sampled ensemble, we first obtain the mass parameter 
mi from C'mix(t) defined previously, and then extract /i 
through a joint fit to Eq. (flUl) with mi fixed. After that, 
we quote the jackknife error of fi as the statistical error. 

For the overlap fermion, the chiral symmetry dictates 
that Zy = Za for the local currents and it has been 
verified in the non-perturbative renormalization . We 
have calculated Za from the Ward identity for a few bare 
quark masses and they are listed in Table m where one 
can see that the quark mass dependence of Za is mild 
and the chiral limit value is Za = 1.1045(8). 


From Zy = Za and the bare decay constant /i, the 
renormalized decay constant fp is obtained at different 
m^, as illustrated in Fig. [H where one can see that fp is 
not linear in m^ throughout the pion mass range. For¬ 
tunately, we have several m,rS very close to the physical 
pion mass, which facilitate us to do an linear interpola¬ 
tion in m^ in the neighborhood of the physical pion mass. 
The final result of fp at the physical point is 

fp± = 208.5 ± 5.5 ± 0.9 MeV, (11) 


where the first error is statistical and the second is the 
combined uncertainty of Zy, the scale parameter a~^, 
and the approximated wall-wall correlation function. 

Experimentally, fp± can be derived from the branching 
ratio of r lepton decaying into p~Vt, 


Br = 


87r 


2 r2 

-rrmpfp- 




( 12 ) 

The branc hing ratio Br(r —> Vr) is measured to be 

25.52(9)% [IJ. After subtracting the 0.30(32)% non- 
p branching fraction, one has the branching fraction 
Br(r —7> P~Vr) = 25.22(33)%. Using this value along with 
the latest PDG2014 values of Gp, |I4d|, rur, Tr, and rup, 
one can get the experimental value = 209.4 ± 1.5 

MeV, where the error is predominantly given by the un¬ 
certainty of the above branching fraction. Our result 
Eq. (fTTl) and are in good agreement. 

To summarize, we use overlap valence fermion for sev¬ 
eral quark masses on the Nf = 2-1-1 domain-wall fermion 
configurations generated by the RBC&UKQCD Collab¬ 
oration. The light sea quark is at the physical point and 
the spatial extension is 5.5 fm. As a benchmark, we first 
calculate the decay constant of pion and obtain a value 
/,r = 131.3(6) MeV at the physical pion mass, which is 
consistent with RBC & UKQCD’s result/,r = 130.2(9) 
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MeV. We extract the spectrum of the 7=1 vector chan¬ 
nel from the Coulomb gauge fixed wall-source correlation 
functions. As such the p mass is precisely determined to 
be rup = 775.9 ± 6.0 ± 1.8 MeV at the physical pion 
mass. We also propose a strategy to reduce the noise of 
the wall-wall correlation functions of hadrons, through 
which the leptonic decay constant of p, fp, is determined 
to be 208.5 ± 5.5 ± 0.9 MeV at the physical which 
agrees well with the value fp = 209.4 ± 1.5 MeV derived 
from the process r —>■ pvT. From the slope of rUp vs. 
we obtain the valence irp sigma term from the Feynman- 
Hellman theorem which gives = 9.82(6) MeV. Fu¬ 
ture results at different lattice spacings are needed for 
the continuum extrapolation. 
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